HEPHY-PUB 812/05 
hep-ph/0510372 
October 2005 



EXACT-PROPAGATOR INSTANTANEOUS 
g ; BETHE-S ALPETER EQUATION FOR 

^ : QUARK-ANTIQUARK BOUND STATES 

-t— > 
o 

o 

; LI Zhi-Feng 

Institut fiir Theoretische Physik, Universitat Wien, 
^ I Boltzmanngasse 5, A-1090 Wien, Austria 

cn 
o 



in 
o 



Wolfgang LUCHA* 

Institut fiir Hochenergiephysik, 
Osterreichische Akademie der Wissenschaften, 
Nikolsdorfergasse 18, A-1050 Wien, Austria 



Franz F. SCHOBERLt 

jr Theoretische Physik, Universi 
Boltzmanngasse 5, A-1090 Wien, Austria 



X 

^ • Institut fiir Theoretische Physik, Universitat Wien, 



Abstract 

Recently an instantaneous approximation to the Bethe-Salpeter formalism for the analysis 
of bound states in quantum field theory has been proposed which retains, in contrast to the 
Salpeter equation, as far as possible the exact propagators of the bound-state constituents, 
extracted nonperturbatively from Dyson-Schwinger equations or lattice gauge theory. The 
implications of this improvement for the solutions of this bound-state equation, that is, the 
spectrum of the mass eigenvalues of its bound states and the corresponding wave functions, 
when considering the quark propagators arising in quantum chromodynamics are explored. 
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1 Introduction 

The more than half a century old Bethe-Salpeter formalism [1] constitutes a relativistically 
covariant framework within the realms of quantum field theory for the description of bound 
states from first principles. The Bethe-Salpeter equation controls a bound-state amplitude 
encoding, in momentum space, the distribution of the relative momenta of the bound-state 
constituents. Within elementary particle physics this formalism has been widely applied to 
quantum electrodynamics (QED) and quantum chromodynamics (QCD). Unfortunately it 
faces problems of interpretation and of our ignorance of the full interaction kernel in QCD. 

Thus, simplifications of the Bethe-Salpeter equation in form of some three-dimensional 
reductions are highly desirable. The most famous among all proposals is known as Salpeter 
equation [2] . Its formulation, however, is based on assuming all bound-state constituents to 
interact instantaneously and to propagate as free particles; the latter circumstance renders 
hard to implement effects such as spontaneous chiral-symmetry breaking, crucial for QCD. 

In view of this, an instantaneous Bethe-Salpeter equation which incorporates the exact 
form of the propagators of the bound-state constituents (to the utmost conceivable extent) 
has been derived recently [3]; this improved bound-state equation reduces, of course, to the 
Salpeter equation, upon approximation of the exact propagators by their free counterparts. 
For any description of hadrons as bound states, the exact quark propagators conforming to 
the QCD Dyson-Schwinger equations are relevant. This work is devoted to the study of the 
consequences of introducing exact quark propagators in this instantaneous Bethe-Salpeter 
equation. The most dramatic effect observed is a significant diminution of the level spacing. 

The paper is organized as follows. Section 2 sketches the derivation of our instantaneous 
Bethe-Salpeter equation with exact propagators of the bound-state constituents presented 
in Ref . [3] . Approximating all interactions entering in the Bethe-Salpeter equation by their 
static forms but retaining, as far as possible, exact propagators yields [3] a generalization of 
Salpeter's equation [2], with momentum-dependent masses of the bound-state constituents 
and with normalization factors of their exact propagators multiplying all interaction terms. 
Section 3 introduces the exact light-quark propagators obtained within QCD as solution of 
the Dyson-Schwinger equations. This infinite tower of coupled integral equations calls for a 
truncation [4-24] which must not be in conflict with the relevant Ward-Takahashi identity. 
Section 4 summarizes the technique developed for finding the solutions of an instantaneous 
Bethe-Salpeter equation by first reducing it to a coupled set of radial equations [25-28] and 
then converting it to a matrix eigenvalue problem [29-32] . Some implications of taking into 
account exact instead of free propagators of bound-state constituents are analyzed in Sec. 5 
by application of the entire formalism to a linear confining interaction. Section 6 scrutinizes 
our findings. Appendix A recalls the Hilbert-space basis required for the matrix conversion. 

2 Instantaneous Bethe-Salpeter equation with exact 
propagators 

The derivation of the exact-propagator instantaneous Bethe-Salpeter equation [3] parallels 
the (three-dimensional) reduction of the Bethe-Salpeter equation [1] to the free-propagator 
Salpeter equation [2] . It may be achieved by several slightly different but equivalent routes. 
In the framework of the Bethe-Salpeter formalism, a bound state |B(P)) of momentum 
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P and mass Mb, composed of a fermion and an antifermion described by the field operators 
■01 (xi), '^2(3:^2), resp., is represented, in momentum space, by the Bethe-Salpeter amphtude 

^(p) =exp(iPX) j d^x exp(ipx) (0|T(V^i(xi) 7/^2(2^2)) |B(P)) . 

Here, X and x denote the center-of-momentum and relative coordinates of the two-particle 
system while P and p label the total and relative momenta of the bound-state constituents. 
This Bethe-Salpeter amphtude ^ has to satisfy the homogeneous Bethe-Salpeter equation 

=(^"^1^1) Jd'qK{p,q)^{q)S2{-p2) • (1) 

The dynamical ingredients of this equation of motion are the exact propagators Si{p) of the 
two bound fermions i = 1,2 (with individual momenta pi, P2) and the interaction kernel K, 
a fully truncated 4-point Green function which encompasses all Bethe-Salpeter irreducible 
Feynman diagrams for two-particle into two-particle scattering and depends on the relative 
momenta of initial and final scattering states, p and q, as well as on the total momentum P. 

The instantaneous approximation to this formalism assumes that the kernel K depends 
just on the spatial components p and q of the relative momenta p and q: K{p, q) = K{p, q) . 
Its application reduces Eq. (1) to the instantaneous version of the Bethe-Salpeter equation 

= ^ [dpo SM I{p) S2{-P2) (2) 
Ztt J 

for the Salpeter amplitude (defined by integration of ^'(p) over the time component po of p) 

Hp) = ^J cipo*(p) ; 

here the term involving the by assumption now instantaneous interaction is abbreviated by 

nP) = -^Jd'qK{p,q)Hq). 

The fermion propagator Si{p) is the solution of the fermion Dyson-Schwinger equation. 
By Lorentz covariance S'j(p) is defined by merely two (Lorentz-scalar) functions Mj(p^) and 
Zj(p^), in QCD referred to as the quark wave-function renormalization and mass functions: 

In the course of the derivation [3] of a generalization of the Salpeter equation towards exact 
propagators of all bound-state constituents, two of the present authors (W. L. and F. F. S.) 
assumed these functions, Mj(p^) and Zi{p^), to depend only on the spatial components p of 
the momentum p. This allows to substitute Miij?) by Mi(p^) and ^i(p^) by Zi{p'^) in Si{p). 

Then the integral in Eq. (2) over the time component po can be easily given analytically. 
Introducing the one-particle energy Ei (p) , the (generalized) Dirac Hamiltonian Hi (p) , and 
the energy projection operators Af (p) for positive or negative energy of particle i = 1, 2 by 

Ei{p) = + ^=1,2, 

H,{p) = 7o[7-p + M,(p2)] , ^ = l,2, 
^ _ Hp)±H,{p) 
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few rather standard manipulations yield [3] our instantaneous Bethe-Salpeter equation for 
fermion-antifermion bound states, with exact propagators of the bound-state constituents: 

m - ZAP.) Up2) 1^ P,-E,{v.)-E,{v,) - Po + Mp,) + Mp,) ) ■ 

Prom this, each amplitude has to satisfy the two constraints A.f{pi) (^2) = 0. 

With little effort our bound-state equation (3) may be rephrased as eigenvalue problem, 

i/l(Pl)$(p)-$(p)i/2(P2) 

+ Z,{pI) Z^ipl) [A+(Pi) 70 Hp) 70 A2 (P2) - Ar(Pi) 70 1{p) 7o ^t{P2)] 
^Po^p) , 

with the bound state's energy Pq or, in its rest frame P2 — — Pi, the mass Mb as eigenvalue. 

The Salpeter equation [2] is obtained by one further step of simplification. Its derivation 
assumes, in addition to the instantaneous approximation for K, that each exact propagator 
in Eq. (1) can be replaced by the propagator So{p, rrii) of a free particle of effective mass mf 

Si{p) = So{p, rrii) = , ■ = i ^ ^ ^ . , i = 1, 2 . 

p — mi + le p^ — mf + ie 

Thus the Salpeter equation is recovered from Eq. (3) in the limit Mi{p'^) rrii, Zi{p'^) 1. 
The exact-propagator instantaneous Bethe-Salpeter equation (3) generalizes [3] Salpeter's 
equation by replacing rrii by Mi{p^) and introducing factors Zi^p^) in the interaction terms. 



3 Quark propagator from Dyson— Schwinger equation 

Within QCD the Dyson-Schwinger equation for the quark propagator involves, besides the 
exact quark propagator, also the exact gluon propagator and the exact quark-gluon vertex. 
The Dyson-Schwinger equations governing the two latter Green functions couple the quark 
Dyson-Schwinger equation to the infinite hierarchy of Dyson-Schwinger equations. Hence, 
a tractable problem can only be defined by some truncation of this set of integral equations. 

For the present investigation we employ the so-called "renormalization-group- improved 
rainbow-ladder truncation" scheme [4-24] applied to the quark Dyson-Schwinger equation 
and the meson Bethe-Salpeter equation. In this specific scheme the exact gluon propagator 
and the exact quark-gluon vertex are replaced by their (perturbative) tree-level forms. The 
truncation is consistent with the preservation of the axial- vector Ward-Takahashi identity; 
this is important for all questions related to chiral symmetry and its dynamical breakdown. 
In this model, all the dynamical information is encoded in some effective coupling strength. 

Viewed as function of the involved momentum transfer squared this effective coupling is 
characterized by two main features. In the ultraviolet region it approaches the perturbative 
behaviour of the strong fine-structure constant, incorporating thereby asymptotic freedom. 
In the infrared region it exhibits the significant enhancement demanded strongly by studies 
of the Dyson-Schwinger equation satisfied by the exact gluon propagator. In the particular 
Ansatz for this effective coupling strength proposed in Ref. [4] this infrared enhancement is 
represented partly by the integrable singularity of a momentum-space S function, partly by 
a finite- width approximation to this S function. This constitutes the "model of our choice." 
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In the comprehensive study presented in Ref . [4] , the exact quark propagators emerging 
from this truncation model are found as numerical solutions of the quark Dyson-Schwingcr 
equation by fitting main properties of the tt- and K-meson system. Like many treatments of 
Dyson-Schwinger equations the analysis of Ref. [4] has been performed in Euclidean space, 
implying that the quark propagators are obtained as Euclidean-space Schwinger functions. 
Within both QED and QCD, the analytic structure of the exact fermion propagators is still 
the subject of intense investigations which did not provide a definitive conclusion until now 
(cf., for instance, Refs. [23,33,34] and references therein). In order to proceed, we must thus 
assume that the necessary analytic continuation from Euclidean to Minkowski space makes 
sense, at least for the quark propagators. In this case, the numerically computed functions 
M(p^) and Z{p^) in the propagator of the hght u- and d-quarks may be represented (rather 
accurately in a range of spacelike momenta) in analytical form by the parametrizations [35] 



b d 
the values of the parameters a, b, c, d, rrio are fixed by interpolation of the numerical results: 

a = 0.745 GeV , b = (0.744 GeV)^ , mo = 0.0055 GeV , 

c = 0.545 , d = (1.85508 GeV)^ . (4) 

These propagator functions read in their "p^ — 0" approximation required for an analytical 
formulation of the bound-state equation (3) proposed [3] as a generalized Salpeter equation 

M(p') = + m„ , Z(p') = l ^-y. (5) 

The behaviour of the quark propagator functions Mij)^) and Z{p^) as functions of is 
depicted in Fig. 1. For light quarks the mass function M{p^) is, of course, dominated by the 
nonperturbative mechanism responsible for dynamical chiral- symmetry breaking. Starting 
at M (0) = 0.7505 GeV, M(p^) drops in the vicinity of = (0.57 GeV)^ by more than two 
orders of magnitude, in order to approach in the limit — > cxo the (comparatively tiny but 
still nonvanishing and hence explicitly chiral-symmetry breaking) current light-quark mass 
mo = 0.0055 GeV. In contrast to such drastic variation, the wave-function renormalization 
function Z{p^) exhibits an only rather moderate dependence on p^. With increasing values 
of p^, Z{p^) rises slowly from ^(0) = 0.455 to its asymptotic value Z{p^) — > 1 for p^ oo. 

Within the "renormalization-group-improved rainbow-ladder truncation" scheme, it is 
by no means mandatory to implement, as done in the model studied in Refs. [4-7,10,11,24], 
the infrared enhancement in the effective interaction by the sum of an integrable 5 function 
singularity and its finite- width approximation. The results of the investigations reported in 
Refs. [8,9,13-18,20,23] demonstrate that propagator functions of very similar shape will be 
obtained in a model in which the infrared enhancement of the effective-interaction coupling 
required by hadron phenomenology is provided only by the finite-width representation [12]. 

Moreover, the predictions for the propagator functions M{p^) and Z{p^) of both models 
[4,8] for the effective coupling in the quark Dyson-Schwinger equation exhibit a remarkable 
qualitative and quantitative agreement with the results produced by lattice gauge theories. 
A recent unquenched lattice calculation of the quark propagator in Landau gauge involving 
two degenerate light (u/ d) and one heavier (s) dynamical quarks may be found in Ref. [36]. 
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Figure 1: Mass function M{p'^) (a) and wave-function renormalization function Z{p'^) (b) 
of the exact propagator of the (light) u- and d-quarks obtained by numerical solution of the 
quark Dyson-Schwinger equation in the "renormalization-group-improved rainbow-ladder 
truncation" model of Ref. [4] as represented by the simple interpolation of Eqs. (5) and (4). 
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4 Pseudoscalar fermion— antifermion bound states 



Now follow the path, paved in Refs. [25-32], of the transformation of bound- state equations 
for Saipeter amplitudes to matrix eigenvalue problems fixing their radial components. 

Consider, as the perhaps simplest example, fermion-antifermion bound states of spin J, 
parity P = (—1)'^''"''^ and charge- conjugation quantum number C = (— l)"^. In spectroscopic 
notation such states are labeled by ^ Jj. Because of the constraints Af (p^) (P2) = 

the general expansion of the Saipeter amplitude over a complete set of Dirac matrices 
involves not the full 16 but only eight independent components. For our Vj states only two 
of the latter, $1 (p) and $2 (p) , are relevant. With our notation for one-particle energy E{p) 
and (generalized) Dirac Hamiltonian H(p) introduced in Sec. 2 the corresponding Saipeter 
amplitude $(p) reads in the center-of- momentum frame of the particle-antiparticle system 



*(p) 
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Without loss of generahty but for definiteness, focus to exactly the same physical system as 
studied in Refs. [29-32] . Aiming at the description of mesons with the quantum numbers of 
the pion discuss quark-ant iquark bound states of spin J = 0, that is, pseudoscalar states of 
spin-parity-charge conjugation assignment J^*" = 0""*". Assume the kernel K{p, q) to be of 
convolution type with Dirac structure that of time-component Lorentz- vector interactions: 
K{p, q) — K{p—q) — V{p—q) 7°(8)7° with V{p—q) any Lorentz-scalar potential function. 

Upon factorizing off all dependence of the Saipeter amplitude on angular variables 
the exact-propagator instantaneous Bethe-Salpeter equation (3) may be reduced [25,27] to 
a set of coupled equations for the radial factors of all the independent Saipeter components. 
For bound states composed of particle and corresponding antiparticle we clearly have, with 
p = IpI, Zi{p^) = Z2(p2) = Z{p^) and Af i(p^) = M2(p^) = M(p2), which will also enter in 
El (p) — E2{p) — E{p) = ^Jp^ -\- M'^{p^). The set of coupled equations governing the radial 
functions $i(p) and $2(1') in our independent Saipeter components and $2(p) of ^Sq 

states, respectively, can be simply derived by, for instance, "dressing" Eq. (1) of Ref. [29] or 
Eq. (1) of Ref. [31] by insertion of the appropriate factors Z(j?) in all interaction terms and 
by replacement of all constant constituent masses m by the relevant mass functions M{p^): 



2E{p) $2(p) + ^'(p' 



(2^ 



Vo{p,q) $2(g) = Mb$i(p) 



2E{p)^^{p) 
+ Z\p') J 



dq q^ 



p 



E{p) 



E{q) E{p) 



Vi{p,q) 



E{q) 



(6) 

$i(g) = Mb$2(p) . 



The configuration- and momentum-space representations of any radial function are related 
by Fourier-Bessel transformations which involve spherical Bessel functions of the first kind 
jn{z) [n = 0, ±1, ±2, . . .) [37] ; as a kind of reminiscence of these, the interaction V{p—q) in 
the kernel K{p—q) enters here in form of some static potential V{r) in configuration space: 

00 

VL{p,q) = ST^jdrr'^jL{pr)jL{qr)V{r), L = 0,1,2,... . 
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The particular structure of the set of equations (6) allows to find its solutions [29-32] by 
inserting one of these relations into the other and obtaining an eigenvalue equation for M|: 



AE\p) ^2{p) + 2Z\p')E{p) 



Aq 

(2^ 



Vo{p, q) $2(g) 



+ 2 



Z^{p^) f dqq 



+ Z 



dq q^ 



M{p')M{q')Vo{p,q)+pqVi{p,q) ^q) 



M{p'^) 
(27r)2 [E{^ 



+ 



p 



M{q^ 

E{q) ' E{p) 



ViiP^l) 



X 



oc 



Vo{q,k) ^,{k) 



(7) 



By expansion over the basis for radial functions summarized in Appendix A we convert this 
integral equation to an equivalent matrix eigenvalue problem which can be diagonalized by 
standard means. As noted in Sec. 2, Salpeter's equation corresponds to the frcc-propagator 
approximation Z{p'^) = 1 and M{p^) = m. Thus the studies reported in Refs. [29-32] could 
get the matrix, for a large class of interactions, in algebraic form. Due to the presence of the 
true quark propagator functions Z(p'^) and M^p"^), in general this is no longer possible here. 
Upon construction of one Salpeter component, $2(p), as solution of Eq. (7), its companion, 
^i{p), follows, for Mb ^ 0, immediately from the first of the two coupled equations (6). For 
vanishing bound-state mass, Mb = 0, Eqs. (6) decouple and are thus solved independently. 



5 Linear confinement: results and discussion 

Let us eventually apply the formalism developed in Sees. 2 through 4 to a confining (static) 
potential of linear shape, V{r) = Xr, with slope A = 0.2 GeV^. For confining interactions a 
time-component Lorcntz-vcctor structure of the kernel appears to be free of all the stability 
problems encountered by solutions found for a kernel of Lorentz-scalar structure [27,38,39]. 

Our first goal is to analyze the effect of the dynamical generation of quark masses on the 
solutions of the instantaneous Bethe-Salpeter equation (3) with exact propagators. To this 
end. Fig. 2 compares, for the three lowest positive-norm J-^*^ = O""*" bound states, the mass 
eigenvalues Mb of Eq- (3) for exact light-quark propagators with mo = (corresponding to 
the chiral limit of QCD) with those of a Salpeter equation for massless constituents [29,30]. 

The chosen bases for the Hilbert space L2{R^) of all with the weight function w{r) — 
square- integrable ( "radial" ) functions on the positive real line introduce one additional 
degree of freedom, by allowing the basis states to depend on a variational parameter > 0. 
As a basis, these vectors constitute a complete orthonormal system for any particular value 
of jji. Therefore, as long as relying exclusively on expansions over the full set of basis vectors 
our results may be expected to be independent of Necessary truncations of expansions to 
a finite number of basis vectors will induce a certain amount of //-dependence of the results. 
However, a reasonable technique involving expansions should exhibit stability with respect 
to the increase of the number of basis vectors. If taking into account a large enough number 
of basis vectors, by reducing the dependence on // some "region of stability" should emerge. 
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Figure 2: Bound-state masses Mb (as functions of our variational parameter /i) of the three 
lowest-lying (positive-norm) J^'-^ = O""*" eigenstates of our exact-propagator instantaneous 
Bethe-Salpeter equation (3) for the propagator parametrization (5) but mo = {full lines) 
and of the Salpeter equation for massless constituents {dashed lines), for a time-component 
Lorentz- vector kernel representing some linear potential V{r) = A r of slope A = 0.2 GeV^. 
These results arise from diagonalization, for given values of //, of 50x50 matrices equivalent 
to our exact-propagator equation but, mimicking Ref. [29] , only 15 x 15 matrices equivalent 
to Salpeter's equation and in both cases 50 terms in all intermediate-step series expansions. 
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With respect to the mass Mb of a given bound state, such a region of stabihty manifests 
itself in form of a plateau where the numerical value predicted for Mb is constant over some 
nonvanishing range of jj,. Beyond doubt, the formation of these plateaus is obvious in Fig. 2. 
Such extrema of Mb disclose the bound states. Compared with the free-propagator results, 
the ground state of Eq. (3) is higher but its radial excitations are lower, to the effect that all 
level spacings are significantly smaller if using exact propagators. In general, it is, of course, 
not possible to compensate for these shifts by some change of the parameter values entering 
in the interaction kernel. For instance, a reduction of the slope A of the linear potential by a 
factor 2 to A = 0.1 GeV^ lowers the ground-state energy eigenvalue of Eq. (3) to the level of 
its Salpeter-equation counterpart but simultaneously diminishes the level spacings further. 

Table 1 lists the masses Mb of the three lowest-lying (positive- norm) J^'^ — bound 
states calculated from the exact-propagator instantaneous Bethe-Salpeter equation (3) for 
the full (mo 7^ 0) parametrization (5) and (4) of the exact light-quark propagators. Because 
of the relative smallness of the (explicitly chiral-symmetry breaking) current mass ttiq these 
eigenvalues are, of course, less than some 0.5% larger than the corresponding "chiral-limit" 
values forming the stabihty plateaus discernible in Fig. 2. These mass values are confronted 
in Table 1 with the corresponding mass eigenvalues Mb of the Salpeter equation, computed 
by assuming, for the constituent mass m of the light u- and d-quarks entering their effective 
propagators, the typical value of m = 0.336 GeV frequently adopted by nonrelativistic and 
relativistic constituent quark models to describe hadrons as bound states of quarks [40,41]. 
Raising in the Salpeter equation the constituent quark mass m from m = to the canonical 
value m — 0.336 GeV shifts the masses Mb of the three lowest states by more than 0.1 GeV 
towards larger values. The net result of this is a further increase of the discrepancy between 
the level spacings predicted by the exact-propagator equation (3) and the ones arising from 
the free-propagator Salpeter equation with some kind of appropriate or reasonable effective 
mass of the constituent quarks. Therefore, we are forced to conclude that any neglect of the 
proper behaviour of the momentum-dependent quark mass, M(p^), by approximating it by 
a constant constituent mass m is, at least for the light u- and d-quarks, rather questionable. 

Table 1: Bound-state masses Mb (in units of GeV) for the three lowest-lying positive-norm 
J^*" = 0^+ eigenstates (denoted by I^Sq, 2^So, 3^So, in usual spectroscopic notation) of our 
exact-propagator instantaneous Bethe-Salpeter equation (3) with light-quark propagators 
parametrized by Eqs. (5), (4) and of the free-propagator Salpeter equation with light-quark 
constituent mass m = 0.336 GeV, for time-component Lorentz- vector kernels representing 
a linear potential V{r) = A r of slope A = 0.2 GeV^, obtained by converting both equations 
to 50x50 matrices and taking into account 50 terms in intermediate-step series expansions. 
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1.813 
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(a) 




Figure 3: Configuration-space radial Salpeter component functions $i(r) (a) and ^2{r) (b) 
for the lowest positive-norm J^*^ = 0""*" bound state of the exact-propagator instantaneous 
Bethe-Salpeter equation (3) with the quark propagator parametrization (5) {full lines) and 
of the Salpeter equation with a light-quark constituent mass m = 0.336 GeV {dashed lines) 
for time-component Lorentz- vector kernels with linear potential V{r) = A r, A = 0.2 GeV^. 
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Figure 3 illustrates for the J^^ — ground state (I^Sq state) of the exact-propagator 
instantaneous Bethe-Salpeter equation (3) with quark-propagator parametrization (5) the 
configuration-space behaviour of the radial Salpeter component functions $i(r) and $2(^)- 
The norm ||$|| of the Salpeter amplitude $(p) for J^*" = 0""*" bound states reads [25,27,29] 

ii^ir - 4 / (If^ ^^^^^ + ^^^^^ ^^^^^^ ' 

this translates for the radial parts $2(p) of the Salpeter components $2(p) to 

ii-^-ir = J dpp^ mip) $2(p) + ■ 

For the plots depicted in Fig. 3, the normahzation has been chosen such that ^2{r) satisfies 

oo oo 

j drr^ |$2(r)|^ = j dpp^ |*2(p)|^ = 1 ; 



the first of Eqs. (6) fixes, after a Fourier-Bessel transformation, the normalization of $i (r) . 
A normalization factor common to $i(p) and $2(p) will then give any desired value to ||$||. 

A closer inspection of the radial Salpeter components $i(r) and ^2{'f) reveals a striking 
similarity to their counterparts found as solutions of the Salpeter equation for a constituent 
quark mass of m = 0.336 GcV. Exact- and free-propagator Salpeter components $i(r) and 
$2(^") show some notable difference only for r < 5 GeV~^ for $i(r) and for r < 2 GeV"^ for 
^2{t). They are hardly distinguishable from each other for r > 5 GeV~^ for both $i(r) and 
^2{t). Thus we feel entitled to expect similar predictions for quantities such as decay rates. 

6 Summary, conclusions, and outlook 

The reduction of the Bcthc-Salpeter formalism to the Salpeter equation requires to assume 
for all bound-state constituents not only an instantaneous interaction but also free-particle 
propagation. Realizing this fact, a bound-state equation that retains the exact propagators 
of the bound-state constituents and generalizes Salpeter's equation has been formulated by 
consequent application of the instantaneous approximation to all exact propagators too [3]. 
Of course, this may be extended to Bethe-Salpeter equations for bound states composed of 
particles that are not, or not all, identical to spin-| fermions as well as to three-dimensional 
reductions [42-45] of the Bethe-Salpeter equation [1] different from Salpeter's equation [2]. 

The present investigation addressed the question of how realistic descriptions of mesons 
as quark-ant iquark bound states within a general instantaneous Bethe-Salpeter formalism 
will be modified by such reinstallment of the exact quark propagators, extracted from QCD 
by analytic continuation from Euclidean to Minkowski space. Interestingly, for the example 
of a specific interaction used already in earlier studies [27-32] , we find a drastic shrinking of 
the level spacings of the bound states while their amplitudes remain practically unchanged. 
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A The generalized Laguerre basis for radial functions 

The Hilbert space L2{R^) of all square-integrable functions on the 3-dimensional Euclidean 
space can be spanned by basis functions each of which is the product of a function of the 
radial variable and of an angular term, the latter being represented by a spherical harmonic 
yimiyi) for the angular momentum £ = 0, 1, 2, . . . and its projection m = —i, —£+1, . . . ,+i 
which all depend on the solid angle fl = {9, 0) and satisfy the orthonormalization condition 



For each value of ^ the radial functions constitute a basis for the Hilbert space L2(-R+) of all 
with the weight function w{r) = square-integrable functions on the positive real line R^. 
The basis functions of L2{R?) in configuration and momentum space are related by Fourier 
transformation. Thus, the configuration-space representation (pf^ (r) and momentum-space 
representation ^ (p) of the radial factors are related by the Fourier-Bessel transformation 



'2 



4\r)=^'J-Jdpp'ji{pr)4\p), t = 0,l,2,..., £ = 0,1,2,..., 





oo 



0f(p) = {-iYJ^Idrr'je{pr)4\r), i = 0,l,2,..., ^ = 0,1,2,.... 



The spherical Bessel functions of the first kind, jniz) {n = 0, ±1, ±2, . . .) [37], are remnants 
of the angular integration. This may be easily deduced with the help of the expansion of the 
plane waves over spherical harmonics in configuration (Q^) and momentum {VLp) space 

oo +1 

exp(ip • a;) = 47r 5^ ^^Upr) yLi^p) yemi^r) ■ 

e=o m=-e 

In terms of orthogonal polynomials of generalized-Laguerre type (for parameter 7) [37], 
L'P\x)-ti-lY(]tl]"4, i = 0,l,2,... , 



which, by construction, are orthonormalized [with weight w{x) = exp(— x)] according to 
Jdxx^exp{-x) Lt'\x) Lf\x) = + ' + S^j , ^, j = 0, 1, 2, . . . , 



our favourite choice of these radial bases is defined in configuration-space representation by 



i2^r+H\ .^^exp(_^^) ^(2^+2)^2 /.r) , z = 0,l,2,... . 



\ T{2£ + i + 3) 

These basis functions involve one positive real variational parameter, with the dimension of 
mass, II. The requirement of their normalizability imposes the constraint fx > 0. Then these 
configuration-space radial basis functions, (j)f\r), satisfy the orthonormalization condition 

00 

j drr^ (j>f\r) cj>f{r) = % , i, j = Q,l,2, . . . . 


Note that the configuration-space representation of our basis functions is chosen to be real. 
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The corresponding momentum-space representation (j)f\p) of our basis functions reads 



(2/x)2^+3i! 



-1) p 



X 



^ (-1)* f i + 2i+2\ T{2e + t + 3) (2/i)* 



^/2£ + t + 3 1 + t ^ 3 p2 



0,1,2,... , 



2 ' 2 ' 2'p2 + ^2y' 
with the hypergeometric series F{u, v; w; z) given in terms of the gamma function F by [37] 



F{u, v; w; z) — 



Tiw) 



E 



V{u + n) V{v + n) z^ 



The momentum-space radial basis functions 0^ "* (p) fulfill the orthonormalization condition 

oo 

I dpp' cl>f\p) <j>f\p) = 5,j , i, J = 0, 1, 2, . . . . 



In momentum space, our basis functions are real for £ = 0, as well as for all even values of i: 
cl>f\p)-4\p) for£ = 0,2,4,... , Vi = 0,l,2,... . 

The virtue of our bases is their analytic availability in configuration and momentum space. 

Mainly for computational convenience, the present investigation makes use of the radial 
basis functions for two values i = and £ = 1 of the angular momentum. Having to deal, in 
momentum-space representation, with the cumbersome hypergeometric series F{u, v] w; z) 
may be avoided by employing simplified expressions equivalent to the above definition [46]: 

I ^1 A i / . , o \ / ^2\-(t+m 

<t>f\p) 



^1 V(-2)*(t + l) I ' + ^ 



i-t 



1 + 



X sin ( (t -I- 2) arctan 



V_ 



lm{(p + i//)2'+3 \p-i(3 + 2i)i4} 



4\p) 



— 1 
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E 



(-2)*/i + 4\ (t + 3)!/i* 



X 



^7^(^ + l)(^ + 2)(^ + 3)(^ + 4) p2 ^ t\ 
— ^ ^ ^ ' sin ( (t + 2) arctan — ] — 

t + 2 y ' 1^1 



i-t J (p2 + ^2-)(t+3)/2 



sin i (t + 3) arctan — 

t + 3 y ' fi^ 



2 ^Jn^ + + 2) {i + 3){i + 4) p^ (p2 + ^,2)3 



X Im|^^-^^^[3/ + 3i(5 + 2i)pV- (5 + 2z)^p//^ -i(5 + 2i)/x^] 
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